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The emerging cosmological picture is of a spatially flat universe composed predominantly of three 
components: ordinary baryons (£Ib ~ 0.05), non-baryonic dark matter (Qjjark ~ 0.22) and dark 
energy (S1a ~ 0.7). We recently proposed that ordinary matter was synthesised from mirror matter, 
motivated by the argument that the observed similarity of Qb and Qoark suggests an underlying 
similarity between the fundamental properties of ordinary and dark matter particles. In this paper 
we generalise the previous analysis by considering a wider class of effective operators that non- 
gravitationally couple the ordinary and mirror sectors. We find that while all considered operators 
imply Qoark = fewxfis, only a subset quantitatively reproduce the observed ratio Qb /Dcart ~ 
0.20. The ~ 1 eV mass scale induced through these operators hints at a connection with neutrino 
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oscillation physics. 
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I. INTRODUCTION 



PACS numbers: 95.35.+d,98.80.Cq,11.30.Er,12.90.+b 
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\Q . A variety of evidence, culminating in the Wilkinson Microwave Anistropy Probe measurements of the cosmic 
£NJ ' microwave background [lj], points to a spatially flat universe composed of approximately 5% baryons (B), 22% non- 
baryonic dark matter (DM) and roughly 70% dark energy. These fractions, while different, are suspiciously similar 
in magnitude. From a theoretical point of view, this is surprising, as one might a priori expect the physics of each 
component to be rather different. This is a new naturalness puzzle. 

The similarity of the positive and negative pressure components, roughly 30% and 70% respectively, has received 
i-S^ ' some attention in the literature. However, the similar magnitudes of Ob and f2_Darfc are also puzzling. (As usual, fix 
denotes the ratio of the energy density of component X and the critical density.) In Ref. |2|, we began an exploration 
O ' of a possible solution. We proposed that the similar baryonic and dark matter densities suggest that the internal 
microphysics of each component are actually similar or identical. Specifically, this will be case if the dark matter is 
■ identified with mirror baryonic matter, B' } The purpose of this new paper is to extend the analysis of Ref. 

A theoretical motivation for mirror matter 0,0 is to retain the full Poincare group, including improper Lorentz 
transformations such as parity inversion and time reversal, as an exact symmetry group of nature, despite the V — A 
character of weak interactions. Mirror matter theories have a G ®G gauge group structure, with ordinary and mirror 
particles assigned to (i?, 1) and (1,-R) representations, respectively. Under parity transformations, a given ordinary 
particle interchanges with its mirror partner. For fermions, a left-handed (right-handed) ordinary particle is paired 
with a right-handed (left-handed) mirror partner. Time-reversal invariance T' follows from the CPT theorem, with 
T'P' = CPT where P' is the exact mirror parity operation defined above. We shall focus on the simplest mirror matter 
model, where G is the standard model gauge group, there are no exotic fermions besides the mirror partners of known 
quarks and leptons, and there is just one Higgs doublet paired with one mirror Higgs doublet. We shall consider only 
the case where the vacuum also respects improper Lorentz transformations, that is, the discrete spacetime symmetries 
are not spontaneously broken |5(. This means that an ordinary particle and its mirror partner have the same mass. 
A large region of Higgs potential parameter space allows this aesthetically appealing scenario [f| . 

Mirror particles interact amongst themselves through their own versions of the strong, weak and electromagnetic 
interactions. The two sectors are inevitably coupled by gravity, but various non-gravitational interactions can also 
connect them. Within the class of renormalisable, gauge-invariant and mirror-symmetry-invariant operators, these in- 
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1 The fig ~ VI Dark issue was very briefly touched upon by Bento and Berezhiani also in the context of mirror matter models. Note, 
though, that their mechanism produces Qb = f° r the symmetric mirror model. 
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teractions are photon-mirror-photon kinetic mixing and Higgs-boson-mirror-Higgs-boson coupling through a <jy (fxj) 1 ' <f>' 
term (we denote mirror fields with a prime) p| . The implications of these interactions have been explored in several 
papers [El El- 

The mechanism under study here employs dimension-5 effective operators of the form lepton-Higgs-lepton' -Higgs' to 
couple the ordinary and mirror sectors. It also requires an assumption about the initial lepton and baryon asymmetries 
(prior to reprocessing). In Ref. Q, a particular such operator with natural initial asymmetries (see below for a review) 
was found to yield 

--5- ~ 0.20 - 0.21, (1) 

* 'Dark 

in excellent agreement with the allowed range suggested by WMAP 1] 0.20 ± 0.02. 

We would now like to know if this encouraging quantitative success is unique to the effective operator we happened 
to consider in Ref. j2j. The purpose of this paper is to examine a wider class of operators. We shall find that while all 
selected operators yield Qb /^Dark < 1, only a subset are quantitatively consistent with the observations, and that 
the operator of Ref. Q is not uniquely favoured (from this cosmological perspective). 

The next section reviews the dynamical framework of Ref. and the section after that presents our calculations. 
We end by discussing our results and drawing conclusions. 



II. REVIEW OF THE DYNAMICAL FRAMEWORK 



A fundamental feature of mirror matter cosmology is that the temperature of the mirror sector, T', is expected to 
be different from that of the ordinary sector, T, at least during certain epochs. For the late epoch during which big 

bang nucleosynthesis takes place, T'/T ~ 0.5 should hold in order to constrain the expansion rate of the universe at 
that time. From a later period still, large scale structure formation with mirror DM suggests a slightly more stringent 
constraint, T'/T ~ 0.2 0]. The smaller this ratio is, the more does mirror DM resemble standard cold DM during 
the linear regime of density perturbation growth (they obviously must differ in the non-linear regime because mirror 
DM is self-interacting, chemically complex and dissipative) jjj. For earlier epochs, however, we have no observational 
constraints that T'/T must satisfy. If £lDark = Q'b-i our fundamental hypothesis, then the inequality of Q! B and Qb 
strongly suggests that the temperatures were also different for at least part of the time during baryogenesis. 

A temperature difference can be created, for example, through inflation Q. Imagine that there is an inflaton 
and a mirror inflaton, and that inflation is seeded by a fluctuation through whichever of these fields the fluctuation 
"favours" . Upon reheating, the macroscopic ordinary-mirror asymmetry generated by the amplified fluctuation will 
translate into T" ^ T provided the two sectors are weakly enough coupled to each other. The subsequent evolution 
of T'/T depends on the precise nature of the ordinary-mirror coupling terms. We emphasise that the macroscopic 
temperature asymmetry is not at all inconsistent with an exactly symmetric microphysics. 

We can now review the main dynamic events our scenario requires: 

Step 1. Suppose that reheating after inflation leaves a universe with 

T' » T, (2) 

that is, the universe is totally dominated initially by mirror matter. This would happen if inflation was driven by the 
vacuum energy of a mirror inflaton. 

Step 2. Around a certain temperature, T' =T\, mirror baryon and/or mirror lepton asymmetries are created. We 
prefer not to specify the mechanism. It might be the out-of-equilibrium decays of the heavy neutral mirror leptons 
as per leptogenesis Q, or something else. No significant ordinary-sector asymmetries are generated - a reasonable 
assumption in the limit of Eq.(J2Jl. 

Step 3. These initial asymmetries are chemically reprocessed. In particular, ordinary asymmetries are created 
through effective dimension-5 operators of the form 

L =(j^) ^L^e'^^' + H.c, (3) 

where £l is a left-handed ordinary lepton doublet, £' R is its mirror partner, (j) is the ordinary Higgs doublet, and 
(j)' is its mirror partner. The indices i,j — 1,2,3 denote the three families. These are the lowest dimension non- 
renormalisable, gauge-invariant operators one can construct out of the fundamental fields. They may be generated 
from renormalisable operators through the exchange of gauge-singlet fermions, hence the notation Mjy. They also 
induce terms m„ = (4>) 2 /Mn in the light neutrino mass matrix. As well as contributing (with sphaleron and other 



3 



effects) to the chemical reprocessing of asymmetries, these interactions also thermally equilibrate the sectors, that is, 
induce T = T' . They operate during the temperature regime, 

M N ~ T ~ T 2 = 10 w {eV/m u ) 2 GeV. (4) 

For T ~ Mff, the effective Lagrangian approach is not valid and the parent fundamental, renormalisable interactions 
are slower than the expansion rate, while for T ~ Ti the effective interactions are themselves slower than the expansion 
rate. The issue explored in this paper is how different family hierarchy assumptions for the (1/Mjv)y affect Qb/Q'b- 
Step 4- A second but relatively brief inflationary episode (or some alternative process) must then occur (beginning at 

some temperature, T3) in order to set up the mild hierarchy, T'/T ~ 0.2 — 0.5, as required for big bang nucleosynthesis 
and the later linear perturbation growth periods. It is interesting that this hierarchy is in the opposite sense to that 
created by the first inflationary episode [see Eq. J5J)]. If the mirror inflaton induces the earlier inflationary burst, then 
it is tempting to ascribe the later burst to a largely failed attempt by the ordinary inflaton to reciprocate. 



III. CALCULATING n B /n Dark 

The main issue is the family structure assumed for Eq. We shall take for simplicity and dcfinitcncss that one 
flavour combination for the effective operators dominates [through having the largest (1/Mjv)y], with all the others 
too small to affect the chemical reprocessing. Given the connection with neutrino mass, we shall assume the relevant 
Mm is of order ((b) 2 /{l eV), since the ~ 1 eV scale is a reasonable upper "bound" on neutrino mass, and is even 
directly suggested by the LSND anomaly [see also Ref. for discussions of ordinary- mirror neutrino mixing] . 

In Ref. Q, we examined only one case, where the i = j = 2 term dominates. Assuming that there is just one eV 
mass term connecting the ordinary and mirror sectors and that it is approximately flavour diagonal, then there are 
just 6 distinct cases to consider (including i = j = 2 to be called case 1 from now on): 

(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

These operators affect chemical reprocessing by constraining the chemical potentials of the species concerned. We 
now need to review how the reprocessing is analysed. 

For temperatures below about 10 12 GeV, QCD 01 an d electroweak non-perturbative processes plus the 
Yukawa interactions for the fermions c, r, b and t are faster than the expansion rate of the universe |l4j. The 
Yukawa interactions for en, ur, cIr, /ir and sr do not become fast enough until the temperature drops below about 
10 10 GeV. Diagonal and off-diagonal weak interactions involving left-handed quarks are also happening rapidly. 
By T 2 ~ 10 10 GeV [see Eq. Q], the selected effective operator, one of Eqs. <|5I1UI1 . has also been inducing rapid 
interactions, affecting the chemical composition and inducing T — T' . At about T — 10 10 GeV ', all the quarks, 
leptons, Higgs bosons and their mirror partners are in thermal equilibrium with distribution functions governed by 
the temperature and chemical potentials for all the involved species. We denote the chemical potential for species 
X by fix (X' by fi x ). The rapid processes listed above relate the /i's. In addition, we impose electric charge or 
hypercharge neutrality, and mirror electric /hypercharge neutrality, for the universe. 
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The chemical constraint equations are 

9/ig + ^""^ nt f — (Electroweak non-perturbative), 

i=i 

3 

6/i 9 — (fi Ui + Hd t ) = (QCD non-perturbative), 
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(r-lepton Yukawa), 



(11) 

plus the corresponding seven equations from the mirror sector. The simplified notation here is: \x q = fj, qiL — ^ q2L — 
fi qaL , where the equalities are enforced by off-diagonal weak interactions; and [ii t = (j,£ iL . 

We are working in the Yukawa-diagonal basis (so u^n becomes the right-handed component of the mass eigenstate 
t-quark after the electroweak phase transition, and so on). 

Then there are one or two more constraint equations, induced by the dimension-5 operator(s): 

(12) 
(13) 
(14) 
(15) 
(16) 
(17) 

So cases 1-3 have 15 constraint equations, while cases 4-6 have 16. With 28 chemical potential variables, this leaves 
13 free /i's for cases 1-3, and 12 for cases 4-6. The free variables correspond to conserved quantities. 

The solution strategy is simply to (i) select the appropriate number of /z's (13 or 12) as independent variables, (ii) 
solve for the remaining chemical potentials in terms of them, (iii) identify the independent, conserved quantities for 
each case and write the associated chemical potential as a linear combination of the independent chemical potentials 
chosen at step (i). Finally, (iv) solve for the baryon/lepton and mirror baryon/lepton asymmetries at T ~ 10 10 GeV 
in terms of the conserved charges. Below 10 10 GeV, other processes come into play in determining the final, low 
temperature values of the baryon and mirror baryon asymmetries via a by now standard procedure (see below) . 

We now examine each case defined by Eqs. Q5I1(J[I separately. In the main text, we shall write down the conserved 
charges and give the expressions for B, L, B' and V at T ~ 10 10 GeV in terms of those charges. In the Appendix, 
we give the algebraic details for case 1 to illustrate the methodology. 
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A. Cases 1 & 2 



Considering first case 1, at T ~ 10 10 GeV there are 13 conserved charges: 

Co = g-B — L2 + 3-B' — L' 2 

Ci = ^B — Li, C\ = 3-B' — L[, 

£2 = 3-B — L3, C 2 = g-B' — L' 3 , 

£3 = L eiR , C 3 = L eiR , 

£4, = L e2R , £4 = L e2R , 

£5 = B U1R - B dlR , C 5 = B' UiR - B' diR , 

£ 6 = B dlR - B d2R , C 6 = B' diR - B' d2R . (18) 

The notation is: Li is the lepton number of family i, B [L — L\ + L2 + L3) is total baryon(lepton) number, B qiR 
is the charge for right-handed quark q from family i, and L eiR is the charge for the right-handed charged lepton of 
family i. The respective quantities from the mirror sector carry a prime. 
In terms of these quantities, 

6 6 

B = «q£q + QLjUj + a'iC'i, 

t=l i=l 
6 6 

t=l i=l 
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TABLE I: The a and coefficients in the expansions of B and L for cases 1 & 2 as defined in Eqs. I19H and 120H . 
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23385 
98276 


n'. - 

"4 — 


15801 
98276 


Pi — 


59567 
98276 


ql _ -20381 
" 4 98276 


Q-5 = 


-963 
24569 
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TABLE II: The a and coefficients in the expansions of B and L for case 3. 



«o = # 




P0 - 158 




Q l = Q 2 = il 


a'i = a' 2 = 


01 = 02= 


Pi - A* - 158 


ft 3 = °4 = 


«3 = q4 = t| 


#3 = 04 = $pf 


03= 04 = Ifil 


flr - =46 

U5 — 869 


rw' 112 
"5 — 869 


a 88 
V 5 — 237 


ai _ -70 
^5 237 


016 = OI5/2 


a 6 = a's/2 


#1 = 05/2 
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Under mirror symmetry, B <-> L <-> L' , d <-> £^ (and £0 -Co)- Hence, 

6 6 

B' = a £o + otjC'i + Q'£j, 

i=l 1=1 
6 6 

£' = foCo + Y.Wi+^fiU (2°) 

i=l i=l 

The values of the a and (3 parameters are given in Table [I] 

Results for case 2 follow from case 1 with the replacements: L<z — ► Li, L' 2 — > in £0 and la (L'jJ — » L2 (L 2 ) in L\ 
(C'±). Because both muon- and electron-Higgs Yukawa interactions are negligible above 10 10 GeV, the a, (3 coefficients 
for case 2 are trivially the same as case 1 which were given in Table [I] 

B. Case 3 

In this case, we again have 13 conserved charges, as in Eq. (|18f) . but with L3 — > L2 (L' 3 — > L 2 ) in £2 (£2)1 an d 
L2 — » L3, L' 2 —> L' 3 in £o- The a and /3 values are listed in Table ITTI 
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TABLE III: The a and f3 coefficients in the expansions of B and L for cases 4 & 5 as defined in Eqs. J2U and (S3). 
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C. Cases 4 & 5 



Considering first case 4, at T ~ 10 GeV, there are 12 conserved charges: 
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In terms of these quantities, 



Under mirror symmetry, B <-> B' , L <-> L' ', £i <-»■ Hence, 

5 



(21) 



5 

i=0 
5 

L = (22) 

i=0 



S' = ^(aid + a'iJCi), 

i=0 

5 

L' = + (23) 



The values of the a and /3 parameters are given in Table ITTT1 

Results for case 5 follow from case 4 with the replacements: L' 2 — ► L[ (Lq, — * L\) in £0 (4) and Li — > L2 
(L[ — > L 2 ) in £1 (£1). Because both muon- and electron-Higgs Yukawa interactions are negligible above 10 10 GeV, 
the a, (3 coefficients for case 5 are trivially the same as those of case 4 (given in Table ^QJ. 



D. Case 6 



This case is similar to case 4, except that L\ — > L3 (L^ — > L' 3 ) in L\ (£[) and L3 — > Li (i 3 — > L^) in £0 (4)- The 
a,/? coefficients are given in Table llVl 
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TABLE IV: The a and (5 coefficients in the expansions of B and L for case 6. 
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IV. RESULTS AND CONCLUSIONS 



The values of the conserved charges, &, CJ i and Co depend on the initial asymmetry generation mechanism. As in 
Ref. 0, we consider, for definiteness, the simple case of non-zero B' and/or L': B' — X , L' = Y" ', B = L = (with 
L' ti = L\ = L'i = Yq/3 and B' UiR = B' d — B' d ). The only nonzero conserved charges are then 

C[ = C' 2 =C = \(X - y ') = Z cases 1 - 3, 

C\ = £' = £ = Z cases 4-6. (24) 
After chemical processing, the baryon and mirror baryon asymmetries at T ~ 10 10 GeV are then 

B = Z(a + a[ + a' 2 ), B' = Z(a + a\ + a 2 ), cases 1 — 3 

B = Z(a + a' + a[), B' = Z(a + a' + a{), cases 4 - 6. (25) 
The ordinary matter/dark matter ratio at that temperature is therefore 
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B 1 
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ba' - 
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(26) 

The value of this ratio changes to some extent at lower temperatures as different chemical processes become 
important. However, at temperatures near that of the electroweak phase transition, T = Tew ~ 200 GeV, the values 
of B and B' depend only on the values of B — L and B' — L' . These charges are separately conserved for T <C 10 10 
GeV, because the interactions in Eqs. I|5I8|) chemically connecting the ordinary and mirror sectors are slower than the 
expansion rate. This yields the well-known relation between B and B — L 15], 

28 , , 
B = -(B-L), (27) 

with an identical relation for B' in terms of B' — L' . Below the electroweak phase transition temperature, there are no 
processes fast enough to further affect B and B' . Thus the final, low temperature value of the ratio B/B' = Hb/^'b 
is simply given (B - L)/(B' - U) evaluated at T ~ 10 10 GeV: 

— - — = — — ~ 0.22, cases 1 & 2; 

\V B a +ai +a 2 — Po — [3i — p 2 256 
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3.25, case 3; 
0.52, cases 4 & 5; 

3.56, case 6. (28) 

If there is some brief period of inflation between T ~ 10 10 GeV and Tew, a s suggested by Step 4, then the results 
depend on when this second period of inflation occurs, T = T3, as well as the subsequent ordinary and mirror sector 
reheating temperatures, Trh and T' RH respectively (with T' RH < Trh required for successful big bang nucleosynthesis 
and large scale structure formation). 

At one extreme, if T3 ~ Tew, then the results of Eq. I|28() hold irrespective of the reheating temperatures. The 
other extreme case, T 3 = 10 10 GeV, allows three outcomes. If T RH > T E w an d T' RH > T E w, then Eq. l|2*S|) again 
holds. In the opposite situation, T R h < Tew an d T' RH < Tew, the final values are simply given by Eq. (|26|l . because 
no further reprocessing can take place. The acceptable intermediate situation, Trh > Tew an d T' RH < Tew, has 
further ordinary sector reprocessing but a frozen mirror sector. For this situation, the final ratios are given by 



tt' B 79 a + ai + a 2 15487 

n B 28 a + a[ + a' 2 - O ~ Pi - P2 _ 182 _ 

fl' B ~ 79 a + ai + a 2 ~ 867 ~ 

n B _ 28 a + a' Q + a[ - (3 - f3' ~ j3{ _ 1578892 



Q' B 79 a + a' + a x 3429785 

n B _ 28ao + a + ai - /3 -Po~P{ _ 440 



0.21, case 3; 

0.46, cases 4 & 5; 



Q' B 79 a + a' Q + a x 813 



0.54, case 6. (29) 



These results all reproduce the qualitative observation that there is more dark matter than ordinary matter. 
Intuitively this is because the FRW universe is born full of mirror matter (under our assumptions), and only some 
of the net mirror baryon/lepton number is chemically reprocessed into an ordinary baryon asymmetry. But it is also 
interesting that there is a subset of effective dimension-5 operators which are quantitatively successful, namely cases 
1-3 [with case 3 marginal unless the circumstances leading to Eq. I|29|l obtain] . Given that the effective operators also 
contribute to the light neutrino mass matrix through ordinary-mirror neutrino mixing, a tentative connection between 
the dark matter problem and neutrino oscillation physics can be made. The connection must be tentative, because 
some important assumptions lie behind our results, especially the microphysical desert between the electroweak scale 
and physics at 10 10-12 GeV (as emphasised in Ref. p|)- 

In conclusion, we have shown that the ratio of baryonic to non-baryonic dark matter, /fldark = 0.20±0.02, inferred 
by WMAPQ, can be quantitatively explained if mirror matter is identified with the non-baryonic dark matter. Our 
explanation involves a set of assumptions about the physics governing the early evolution of the Universe, which 
are not unique but are nevertheless plausible. Our approach also has important implications for neutrino physics, 
suggesting eV scale neutrino masses, which can be tested/constrained from upcoming neutrino experiments such as 
miniBooNE. 



Acknowledgments 



This work was supported by the Australian Research Council. 



APPENDIX A: SOLVING THE CASE 1 EQUATIONS 



We show the algebraic technique for solving case 1. The other cases follow similarly. 

The aim is to use Eqs. <|11[) and (|12fl to extract the total baryon number B, total lepton number L, and their mirror 
matter analogues B' and L' at T ~ 10 10 GeV. These 15 equations reduce the 28 variables, fi q , (j, Ui , /j,^, Hu, fX(j> 
and their primed counterparts to 13 independent variables. The number of independent variables corresponds to the 
number of conserved charges, Eq. I|18(l . The problem at hand is to find B, L, B' and L' in terms of the conserved 
charges. 
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One systematic way of doing this is the following. First identify 13 independent variables. One possible choice is 
the following: fi q , /^ 2 , fj, i3 , Mei, Me 2 , Mdn Md 2 , M^ M* 3 ) Me^ M e2 , and Md 2 - Then use Ec l s - d and 03> to writc 
the 15 dependent m variables in terms of the chosen independent variables. Doing this we have: 





— -9/i 9 - m^ 2 - M«3 










M<A 


= g (-21/A? + 3/idi 


+ 3Md 2 


+ M ei 


+ Me 2 


+ M<? 3 ) 


Md 3 


= « ( 27 Mg ~ 3Md! - 

D 


3Md 2 - 


Me! - 


Me 2 - 


M<? 3 ) 


Me 3 






Mei — 


Me 2 + 


°M« 3 J 


M«3 


= Mu 2 = g (-15/ig " 


1- 3Mdj - 


I- oMd 2 


+ Mei 


+ Me 2 


Mill 


^ /on n 

= g (39/i 9 - - 


9Md 2 - 


Mej - 


Me 2 — 


M<? 3 ) 




— M<" 2 + M</> — M^> 












= ~9Mg ~ Mf 2 - iH 


+ M0 ~ 


M£ 3 






m'«j 


= g ( _21 /4 + 3/4 


+ 3Md 2 


/ 

+ M ei 


+ Me 2 


+ /4) 


Md 3 


= £ ( 27 /4 - 3M dl ~ 


3 Md 2 ~ 


/4 " 


Me 2 - 


/4) 


/4 


- g(21/i;-3/x dl - 


3Md 2 - 




Me 2 + 


5 /4) 


/4 


, 1 , , 
= M« 2 = g (-15M, - 


\- 3M dl - 


+■ 3 /4 


+ M' ei 


+ Ml 2 


M« x 


= i (39^; - 9M dl - 


9/4 - 


a4- 


Me 2 - 


/4) . 



(Al) 

where it is understood that in the m^ 2 equations, m</> an d M^ are to be substituted with the respective righthand 
side's above. In terms of the chosen independent Mi, the baryon and lepton number are: 

3 

B = 6f.i q + (m«, + MriJ 

i=l 

= 12Mg 

3 

L = >^ (2M£, + Me; ) 
i=l 

—29 1 1 5 5 5 

= -^-Mg - 2 /idl ~ 2 Md2 + 6 Mei + 6 Me2 + 6^ 3 ^ A2 ^ 

The next step in the calculation is to write the conserved charges, Ci in Eq. (JTSJ, in terms of the 13 independent 
variables. For example: 2 

A = \ B - L i 



= 2 Mg + \ y^(Mui + MdJ - 2 M£i - Mei 

= 22fi q + 2m£ 2 + 2m£ 3 - Mei , (A3) 
where Eq. IjAlJ) has been used in the last step. The results for the other Ci are given below: 

4 1 1 ...1.2,1, 

Co = llMg - 4m« 2 - -Me 2 - 3 Mei - gM£ 3 - Mdi - Md 2 - 3m 9 + Mdi + Md 2 + gM ei - gM e2 + ^e 3 



2 The conventional definition of, say, the "baryon number of the universe" is the ratio tlb/s where is the net baryon number per unit 
volume, while s is entropy density. There is proportionality factor relating this definition of baryon number to the simple one convenient 
for our application (just the appropriate linear combination of chemical potentials). 
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£3 = Mei 
£4 = fi e . 2 



£5 = - (39/i g - 15/i dl - 9/Z d2 - ^ei - Me 2 - W3) 



1 

6 

£■6 = Mrfi _ Mrf 2 

4 = 29/^ ~ + + Md2 + + ^Me 2 + 2/if 2 + ^3 - Md! ~ Md 2 ~ ~ ^Me 2 + ^/4 3 

£3 = < 

A = a4 

4 = g ( 39 Mg - 15^d! - 9/4 2 - Mei - Me 2 - /4) 

4 = *4 - Md 2 (A4) 

For the given set of values for the conserved quantities, £i and £[, the baryon number can be found by solving the 
identities: 

6 

B = 12/ig = 12 x fi q + x /i£ 2 + x fi ta + ... ee a £ + (a i £ l + a-£-) . (A5) 

i=l 

where it is understood that the £'s are also functions of the 13 independent fi variables [using Ea. (jA4|) ]. By equating 
coefficients of each of the 13 independent variables a set of 13 simultaneous equations for the en's results: 

1 13 

12 = lla + 22ai + -a 2 + —a 5 - 7a[ 

= 2ai - 4a + 2«i 

17 1 1 1 , 

= 2ai - — a 2 - -a s - -an + -a x 
6 6 3 3 

1 111, 

= -ai + -a 2 + a 3 - -a 5 - -a + -a 1 
6 6 3 3 

1 14 1, 

= -a 2 + a 4 - 7«5 - o a o + o a i 
6 6 3 3 

1 5 

= -a 2 - + a 6 - a Q + a x 

1 3 / 

= -a 2 - -a s - a 6 - a Q + a x 

= _3ao + 29ai + iQ4 + yaJ 

1 5 , 17 , 1 , 

= -ao + - y a 2 - -a 5 

1 4,1, , 1 , 

= —an a, H — a 9 + ac 

3 3 1 6 2 3 6 5 

-2 1,1, , 1 , 

= — a ~ g"i + g" 2 + "4 ~ g« 5 

1 5 

= a - a[ + -a 2 - -a' 5 + a' 6 
1 3 

= a - ai + -a 2 - ^a's ~ a '& ( A6 ) 
These 13 equations can easily be solved for the 13 a's; the results are as displayed in tabled 
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A similar procedure with L instead of B in Eq. (|A5|) yields results for the f3's. 
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